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Abstract. - The effective mass of the quasiparticle excitations in quasi two-dimensionai 
systems is calculated analytically. It is shown that the effective mass increases sharply when 
the density approaches the critical one of metal-insulator transition. This suggests a Mott-type 
of transition rather than an Anderson-like transition. The experimental measurements can be 
reproduced in this way without any additional parameter. 



The explanation of the metal to insulator transition (MIT) at low temperatures in quasi 
two-dimensional systems is still a strongly debated task. A critical discussion of different 
approaches can be found in Ref. [1]. The generic feature of MIT transition is the rapid 
change from insulating to conducting behavior when the density is increased very slightly 
at low temperatures. This density driven MIT transitions are usually referred to as Mott 
transitions. The characteristic feature of the Mott-Hubbard transition is that the increasing 
effective mass is the reason for increasing resistivity p = m/e 2 nT while the Anderson scenario 
would assume a vanishing relaxation time nr. It is obvious that this feature characterizes the 
transition rather than the nature of the insulating state itself, see for details [2] . 

In a recent experiment [3] it was shown that the effective mass is increasing sharply when 
approaching the critical density. This would underline the Mott picture rather than the 
Anderson transition. 

Here in this letter we want to substantiate this picture by a quantitative explanation of the 
experimental values of the effective mass. To this end we will introduce a new approximation 
which is based on the large mass difference between transport electrons and scattering impurity 
of donor ions. Our model consists in electrons scattering with heavy ions within the quasi two- 
dimensional gas. In this way we will describe the transition due to Coulomb correlation and 
not the nature of the insulating state itself. Assuming the motion restricted to the x — y plane, 



the Coulomb potential in this cylindrical Fermi surface is V a b(q x , %) = 27re a ef,?i/ J q 2 + q 2 . 



We want to determine the quasiparticle mass which will be compared to the experimental 
results. To this end we will use the standard quasiparticle picture based on the Green's 
function method. As generally known, within this approach the quasiparticle energy and the 
mass of the model are determined by the real part of the selfenergy according to the following 
formulas. 
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Fig. 1 - The selfenergy in terms of the screened potential V = V + VLTV (wavy line). 

First, the quasiparticle energy hui — e k is given as a solution of [E = ReE 71 ] 

k 2 

u E(fc,w)=0. (1) 

2m v ' w 



Second, from Q the effective mass m* follows as 



k_ = de k = k | gE de k | dE = A + g| 
m* 9fc m <9?iw <% 9/c 1 — as 



fe I as 

" " m dk fn\ 



where the arguments have to be put on-shell hto — e k after performing the derivatives. 

In order to proceed we have to know the selfenergy E^. For the scattering of electrons 
on charged ions or holes it is important to consider the particle-hole fluctuations up to any 
order. This is considered in dynamically screened approximation (GW) which expresses the 
selfenergy by a sum of all ring diagrams, given in figure ^ Due to the big mass difference 
between electrons and ions the vertex corrections are suppressed. The experimental data we 
want to compare with are characterized by a relatively low density of electrons such that the 
electron-electron contribution can be neglected. 

To avoid technical complications with the Matsubara frequencies and the subsequent ana- 
lytical continuation to the real time, we will use directly the method of the real time nonequi- 
librium Green's function even for the equilibrium though its strength lies in the description 
of noncquilibrium situations. Within this approach the retarded selfenergy of a species noted 
by the subscript a can be written as 

E?(g,w) = [^-^-r a ( q ,u')-*-r a ( q ,u;) (3) 

J 2n uj — uj' 2 

with the imaginary part 

r a ( 9lW ) = S>(g, W )+E<(g, W ). (4) 

> 

The correlation parts of the selfenergy, E < , in the dynamical screened approximation read 

S<(fc,w) = J^0^^V<(q,u/)G>(k-q,<J-u). (5) 
The screened potential V is expressed via the density fluctuation L according to figure as 



V<{ q ,u) = Y,V{q)l b L<{ q ^) = Y,n<l)lM^ (6) 



b b 

where £ R (q,uj) = 1 — Ylb Vbb(q)^l R {q,uj) is the dielectric function and the polarization or free 
density fluctuation is given by 

n<(«,u;) = Jj0^^G<fa<S)G>(p+q,u/- U ). (7) 
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Since we are on the level of the quasiparticle picture, we can replace in the above for- 

> 

mulas all electron and hole correlation functions G < by their quasiparticle approximation in 
equilibrium 

G<(k,u>) = f<5{e k -hw) (8) 

where — f™, fa — ^ ~ /<f D are given by the Fermi function. 

As the result we get for the selfenergy iJJj} in the dynamical screened quasiparticle approx- 
imation the formula 

> 

Since we want to consider here the screening due to the scattering with heavy impurities 

(or heavy ionic traps or heavy Fermions) we can further approximate this selfenergy. For this 

> 

purpose we rewrite the free density fluctuation or polarization function II < in a different form. 

> 

In terms of the Bose function g(x) = l/(exp (x/T) — 1) the polarization function n < reads 

n£(q,u) = ±2g b {±fiuj)lmU b (q,uj) (10) 

where the imaginary part is 

/dp 
( 2?r fi)2 (f b * " fb,p+q)5 (Cj, ~ € p+q ~hu>). (11) 

A final simplification of the selfenergy can be achieved in the limit of heavy masses of the 
scattering impurities. We are allowed to neglect quantum fluctuations of these impurities 
which are expressed by g(±huj) in i|l(J|) and replace g(±hui) — ► ±T/Tilu. The selfenergy iJSJ 
can then be written as [V 2 b = K a Vf,b] 

£<(fc,c) = 2 / -J% Vaa (q)f< k _ q f ^-5{e k - q -hLJ-hu')Tm ' 



= I (2w Ka(9)Re i 1 ~ em) f> ^-^ - M (12) 

and the real part of the selfenergy becomes 

In order to evaluate the effective mass according to |J2J, it remains to get an expression for 
the dielectric function. To this end we use the zero temperature expansion of the dielectric 
function in quasi two-dimensions [4,5] since the leading temperature dependence is already in 
front of the integral. This means we consider the electrons as degenerated but the impurity 
particle-hole fluctuation as classical here. Changing the integration variables z — q/2pf, we 
then obtain for the real part ifl^l of the dynamical screened selfenergy [x = k/2pf, x = 
huj/4e f ] 

OO 27T OO 2-7T 

e 2 a Tm a f dz f dxj) e 2 a Tm a f f d<f> 



2hp f x J ^ + 1 J ^-l + 2zcos(/)-z 2 2hp f x J J f£-l + 2zcos 



n 
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K 



4x 





-K 



< x a < x 2 

Xq = X 2 

Xo > x 2 

Xq < 



(14) 



with JC the complete elliptic integral of first kind and n a = hx a /2pf. The inverse screening 
length is given by \ a = 2e 2 m a /h 2 . The approximate sign concerns the limit of large n a which 
is justified for the parameters used here since for typical densities of 7 x 10 10 cm -2 we have 

284.9 m* 



y/na/7 x 10 10 cm- 2 m 
In the same way we can evaluate the imaginary part of the selfenergy |0J) with the result 



(15) 



2TT 2 e 2 a Tm a 



K 



(a; + y/xa) 7 



(16) 



The linear temperature dependence appears here from the neglect of the quantum fluctuations 
in the ionic particle-hole fluctuation and is due to screening and should not be confused with 
the standard non-Fermi liquid behavior in the literature. The latter one is more visible in 
the divergency at the Fermi energy xq — x 2 in figure |2 Later we will give also the statically 
screened result which can be considered as a Born approximation of an electron scattering 
with a neutral impurity in the form of a Debye potential. This selfenergy has a finite zero 
temperature limit but the imaginary part does not vanish at the Fermi energy though it is not 
diverging. For a Fermi liquid we would expect a vanishing imaginary part of the selfenergy at 
the Fermi energy. 

The real part E is related to this imaginary T by the Hubert transform The more 
astonishing is the fact that we find here an additional relation 



E(fc,w) = -T(k, 



2mh 



(17) 



which is only valid for this specific type of selfenergy besides the Kramers-Kronig relation 
© . From the plot in figure [21 we see that the excitation are only possible for positive fre- 
quencies since we have calculated the first temperature correction. One can see that this first 
temperature correction has a highly nontrivial frequency behavior far from being Fermi-liquid 
like. 

We can now evaluate the effective mass and compare its dependence on the density with 
the experimental facts in figure|3J Since the real part of the selfenergy (|14|) has the dependence 
on momentum and frequency as E oc l/kF[2muj/k 2 } we see that @ takes exactly the form 



1 ~ ^a(Pf^f) 
The needed on-shell value of the self energy l|14(l simplifies 



m 
m* 



(18) 



2 r 

= 7T — K a 



(19) 
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Fig. 2 
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Fig. 3 



Fig. 2 - The real and imaginary part of the first temperature correction of the selfenergy (1161 and 
11411 versus scaled frequency Xq/x 2 = 2mTiui/k 2 . 

Fig. 3 - The effective mass for temperatures 0.1,0.4, 0.8K from left to right according to 11811 versus 
density. The open circles are experimental values [3] at T = 0.1K. 



with Kn — tt^- 



_ r lXa _ 

2p f 



h Pf 



According to the formulae (|18fl we obtain a sharp increase of the 
effective mass for lower densities (figure |3J) still above the critical one. This is in very good 
agreement with the measurements [3]. 

In figure we compare the formula (|18|l with (|19(l with the experimental values of [3] . 
The agreement is astonishing precise considering the fact that the only parameter entering 
the formulas are the experimental density and temperature. 

The temperature dependence of the effective mass turns out to be linear over a wide range 
since from <|18|) and (|19|l we obtain for small corrections £ 



1 If 1 „ k 2 , 

— « - i-^-s Q fc,— ) 

m* m \ Zef 2m 



k= Pj 



1 / 2 T 

= — 1 - 7T — K a 

m \ e/ 



(20) 



while for large S the equation l |18|l has to be used. This implies a nearly linear temperature 
dependence for the conductivity, a — ne 2 T/m*, if the relaxation time behaves accordingly. 
This has been shown indeed in [5]. The linear temperature dependence has been repeatedly 
reported in the literature both from experimental and theoretical point of view. Numerical 
calculations of Coulomb scattering rates on impurities predict a linear temperature dependence 
of the mobility in silicon inversion layers [6,7]. This was attributed to the collisional level 
broadening in the screening function. Related results have been obtained in Ref. [8] where a 
substantial suppression of the temperature dependence in the screening function was found. 
An analytical investigation of screening in quasi two-dimensional systems was given in Ref. [9] 
where a linear temperature term to the conductivity was found which was confirmed in [5]. 

When the density decrease towards the critical value, the electrons become trapped by 
the charged impurities forming neutral impurities. This is the mechanism of Mott transition. 
A corresponding mass action law will then determine the neutral impurity concentration n,. 
Above the critical density these neutrals are negligible. We suggest here that the experimental 
values in figure are well above this value. 
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For the reason of completeness we give now the contribution to the effective mass if there is 
a scattering with the neutral impurities. We can use a Debye form of the scattering potential 
which results formally from the static approximation in 10 

V ,( a ) - V{q) - 2lT9abh (21) 

if we introduce the scattering strength g a b — e a e-b and the range of interaction ro = 1/k. 
In the second Born approximation [10] the relation between the scattering strength and the 
scattering length reads 

2m 2 / m \ 
«o « --^-gabr 11 + j^gabro I . (22) 

From @ we obtain now 

£<(fc,u;) = 2™, | -J^ Vl (q) 2 f<_ q S(hcj-e^ q ) (23) 
and instead of (|14f) the real part of the self energy reads 



E^) = / *» — (24) 



2- 



The on-shell value xo = £ 2 can be performed analytically to 

k 2 , q 2 ,m n x 2 rii ( 1 Kjarccos 



2m ; " 2^ 2 n a \\x^-K 2 (4x 4 - k 2 ) 3 / 2 

with ^ = h/2ropf. The value at the Fermi energy is 

2 

E *(P/'2^r) A nt f mg ab r \ 2 

— = -4— fc2 (hirtj-1). (26) 

Together with the electron - charged impurities contribution (|20|l we obtain finally 

(ln/ti-1)^. (27) 

When the neutral impurities become significant near the Mott-transition we have to use a 
mass action law to determine Uj. This can be found from a simplified semiconductor model. 
Assuming only the scattering with donor levels, the total donor concentration consists of 
neutral and charged impurities, n = rii + n~l . The neutral donors correspond to electrons 
trapped at the donor level £x>. With the Fermi energy Ep, from the thermal population 
follows [11] 
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where the factor 1/2 comes from the two possible states at each impurity site. The electron 
density in the band is given by the effective conduction band energy E c (n,T), [12] 

m n T %-«cK,t) 

■e T . (29) 



n 



+ 



nh 2 



where E c (n a , T) describes effective conduction band level which becomes density and temper- 
ature dependent by the correlation effects of the electrons. Eliminating the Fermi energy in 
(|28|l and (|29|l one obtains the mass action law [12] 



r e (30) 

n a m a T 

The density dependence of the selfenergy S which determines the effective conduction band 
E c leads to a nonlinear density dependence of this ratio. In principle also the donor level En is 
density dependent. But, we can safely condense both effects into an effective density dependent 
binding energy Eb(n a ,T) — E c (n a ,T) — Ed- With increasing density of the electrons more 
collisions with donors happens and the formation of bound states is favored until a critical 
density is reached where pressure ionization happens. At this critical density the trapped 
states are resolved called here the Mott-transition. The particle distance becomes smaller 
than the Bohr radius. 

As conclusion, the found effective mass here which can describe the experimental values 
seems to support the Mott-transition picture rather than the Anderson scheme, the same 
conclusion is obtained from the calculation of the conductivity [5] where the reader is kindly 
referred to for a more detailed discussion. It should be pointed out however, that the nature of 
the insulating state itself is not clarified by the above consideration though we can describe the 
sharp increase of effective mass by Coulomb correlations. The suggested trapping mechanism 
of electrons on charged impurities [5] is not the only possible one. An alternative is the 
formation of three-particle clusters [13] which is more favourable in high magnetic fields. The 
strong magnetic field dependence seems to underline such idea. Further insight will be gained 
from the calculation of magnetic field dependence than the increase of effective mass which is 
shown merely due to Coulomb correlations. 

The enlightening discussions with Peter Fulde, Pavel Lipavsky, Enver Nakhmedov, De- 
banand Sa, Michael Schreiber and Vaclav Spicka are gratefully acknowledged. 
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